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We discuss the prospects for "tomography" of the intergalactic medium (IGM) at high redshifts using the 21 cm 
transition of neutral hydrogen. Existing observational constraints on the epoch of reionization imply a complex 
ionization history that may require multiple generations of sources. The 21 cm transition provides a unique tool 
to probe this era in detail, because it does not suffer from saturation effects, retains full redshift information, 
and directly probes the IGM gas. Observations in the redshifted 21cm line will allow one to study the history 
and morphology of reionization in detail. Depending on the characteristics of the first sources, they may also 
allow us to probe the era before reionization, when the first structures and luminous sources were forming. The 
construction of high signal-to-noise ratio maps on arcminute scales will require approximately one square kilometer 
of collecting area. 



1. Introduction 

One of the major goals of modern cosmology 
is to understand how the first luminous sources 
in the universe formed and grew into the galax- 
ies and galaxy clusters that we see today. As 
technology has improved in the past four decades, 
the frontiers of our knowledge have extended far- 
ther and farther into the past. We now sit at the 
cusp of exploring the era in which the first gen- 
erations of sources were just beginning to form: 
the end of the cosmological "dark ages." Ob- 
servations of this epoch will teach us about the 
transformation from the smooth, relatively sim- 
ple universe that we see in the cosmic microwave 
background (CMB) to the complex universe that 
we see nearby. 

The hallmark of the era of first objects is 
the "reionization" of the intergalactic medium 
(IGM), when these first luminous sources ionized 
the gas around them and rendered it transpar- 
ent to ultra-violet photons that we observe today 
as optical or infrared light (see [2] and references 
therein). Fortunately, radio photons, including 
those from the CMB, readily propagate through 



both the neutral and ionized IGM at this time, 
and these will help us to measure an enormous 
amount of information about the process of reion- 
ization, the properties of the first sources, and the 
properties of the IGM 50 6.24, . 

Existing techniques have provided tantalizing 
constraints on the epoch of reionization (FOR). 
The most straightforward method is to observe 
the high redshift part of the "Lya forest:" re- 
gions with relatively large HI densities appear as 
absorption troughs in quasar spectra, which pre- 
sumably deepen and come to dominate the spec- 
tra as we approach the reionization epoch. In- 
deed, spectra of z ~ 6 quasars selected from the 
Sloan Digital Sky Survey (SDSS) show at least 
one extended region of zero transmission in- 
dicating that the ionizing background is rising 
rapidly at this time ^21 (but see 03). However, 
the mean optical depth of the IGM to Lya absorp- 
tion is TLya ~ 6.45 X 10^xh[{1 + z)/10]3/2 gg, 
where xh is the global neutral fraction. Even 
a modest xh > 10~^ will therefore render the 
absorption trough completely black; quasar ab- 
sorption spectra can clearly probe only the latest 
stages of reionization. 
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A second constraint comes from the effects 
of the ionized gas on the CMB. The free elec- 
trons Thomson scatter CMB photons, washing 
out the intrinsic anisotropics but generating a 
faint polarization signal. The total scattering 
optical depth t^s is proportional to the column 
density of ionized hydrogen, so it provides an 
integral constraint on the reionization history. 
Recently, the Wilkinson Microwave Anisotropy 
Probe (WMAP) used the polarization signal to 
measure a large Tes, indicating that reionization 
began at 2 > 14 |28I46| . A third (possibly con- 
ficting) constraint comes from measurements of 
the temperature of the Lya forest a,t z ^ 2-4, 
which, taken at face value, require reheating of 
the IGM at z < 10. Since heating is accom- 
plished through photoionization, this would sug- 
gest an order unity change in the ionized fraction 
at z < 10 |47l2()j . though with large uncertainties. 

Taken together, these three sets of observations 
imply a complex reionization history extending 
over a large redshift interval {Az ~ 10). This is 
inconsistent with a "generic" picture of fast reion- 
ization (e.g., [2], and references therein). The 
results seem to indicate strong evolution in the 
sources responsible for reionization, and a de- 
tailed measurement of the reionization history 
would contain a rich set of information about 
early structure formation 44 50 6 24 . In fact, it 
is quite plausible that the ionization history is 
complex, with long periods of "stalling" or even 
two distinct epochs of reionization. Figure ^ il- 
lustrates why. Here we consider the possibility of 
two generations of ionizing sources . In the 
bottom panel we show the ratio of the recombi- 
nation time (assuming the mean cosmic density 
and T = 10* K, typical of ionized gas) to the age 
of the universe (solid curve) . Note that the ratio 
is near unity at z '--^ 10 and decreases at earlier 
times, implying that a significant fraction of the 
gas can recombine if it is ionized by this point. 
The dotted curve shows how xh could evolve 
and roughly satisfy the three observational con- 
straints. The first generation of sources (for ex- 
ample, very massive metal-free stars) ionize most 
of the gas. However, these sources quickly cease 
forming (perhaps because they enrich their envi- 
rons with heavy elements) and most of the gas 



recombines. Full reionization must await the sec- 
ond generation of sources (in our example, Popu- 
lation II stars). We see from the Figure that the 
21 cm structures could be long-lived and complex: 
for example, gas above the mean density may re- 
combine between the two generations while voids 
remain ionized. The top panel illustrates the tem- 
perature evolution of the IGM in this scenario; 
note that reionization corresponds to a sudden 
boost in the temperature. This panel will be dis- 
cussed in more detail below. 
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Figure 1. A cartoon illustrating a double reion- 
ization history. The bottom panel shows the ratio 
of the recombination time and the age of the uni- 
verse (solid curve) as well as xh (dotted curve). 
The top panel shows the kinetic temperature of 
the IGM gas (solid lines), the CMB tempera- 
ture (short-dashed line), and the spin tempera- 
ture (long-dashed line, note that for z < 23 we 
assume this traces Tk)- The thick solid line ne- 
glect all heating mechanisms except for photoion- 
ization, while the thin line shows what could hap- 
pen if X-ray and shock heating were included (see 
§3) . The shaded regions mark the two generations 
of ionizing sources. 
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While such cartoon views are certainly plau- 
sible, detailed models rely on a number of un- 
certain parameters that we would like to mea- 
sure. The optimal reionization experiment would: 

(1) be sensitive to order unity changes in xh (to 
probe the crucial middle stages of reionization), 

(2) provide measurements that are well-localized 
along the line of sight (rather than a single inte- 
gral constraint), and (3) not require the presence 
of bright background sources, which may be rare 
at high redshifts. The most promising candidate 
proposed to date is the 21 cm hyperfine transition 
of neutral hydrogen in the IGM |14ll5j , which ful- 
fills all three of these criteria. 

So long as the excitation temperature Tg of the 
21 cm transition in a region of the IGM differs 
from the CMB temperature, that region will ap- 
pear in either emission (if T5 > Tcmb) or ab- 
sorption (if Ts < Tcmb) when viewed against 
the CMB. Variations in the density of neutral 
hydrogen (due either to large-scale structure or 
to HII regions) would appear as fluctuations in 
the sky brightness of this transition. Because it 
is a line transition, the fluctuations can also be 
well-localized in redshift space. Thus, in prin- 
ciple, high resolution observations of the 21 cm 
transition in both frequency and angle can pro- 
vide a three-dimensional map of reionization. To- 
gether with radio absorption spectra of bright 
background sources (see the contribution by Car- 
illi, this volume), these observations promise to 
shed light both on the early growth of structure 
and on reionization. 

In the following, we will distinguish two ob- 
servational goals. The ultimate goal is to make 
three-dimensional, high signal-to-noise maps of 
the sky in the 21 cm transition. In addition to 
precisely mapping reionization itself, this would, 
for example, allow detailed comparisons between 
the distribution of neutral gas and the ionizing 
sources. Unfortunately, we will find below that 
the expected brightness of the IGM is extremely 
small both in absolute terms and when compared 
to the relevant foreground contaminants. It is 
therefore considerably easier to make statistical 
measurements of the distribution of neutral gas. 
As we have seen with CMB measurements, these 
statistical constraints can be quite powerful. 



In §2 we will briefly review the important 
physics of the 21 cm transition. In §3 we will 
outline the information we can expect to learn 
from 21 cm observations, and in §4 we will de- 
scribe the instrumental requirements for making 
maps and statistical measurements. Finally, in 
§5 we conclude and discuss 21 cm observations in 
relation to other instruments that will be built in 
the next decade. 

2. The 21 cm Transition 

We define 5T{v) to be the observed brightness 
temperature increment between a patch of the 
IGM, at a frequency v corresponding to a redshift 
1 + z — vi^/v (where vq = 1420.4 MHz is the rest 
frequency of the 21 cm transition), and the CMB. 
Assuming the WMAP cosmological parameters 
and small 21 cm optical depth T21 « 1, this 
quantity is |42I33| : 

^T{v) « — T21 (1) 

1 + z 

« 2Z{1 + 5)xhT{-^\ mK, 

where T — {Ts — Tcub)/Ts, Ts is the spin tem- 
perature of the IGM, Xh is the neutral fraction, 
and 6 is the overdensity. All of these local quanti- 
ties should be averaged over the volume sampled 
by the telescope beam (note that we will use xh 
for the global, mass-averaged neutral fraction). 
For reference, one arcminute subtends « 1.9h~^ 
comoving Mpc at z = 10 in the WMAP cosmol- 
ogy, while a bandwidth of 0.1 MHz corresponds 
to w 1.7 Mpc at the same redshift. 

The observability of the IGM clearly depends 
on the spin temperature, which may be written 

mm-- 

j,^ _ TCMB + VcTk + VLyaThya 
1 + yc + ULya 

The second term describes coUisional excitation of 
the hyperfine transition, which couples Ts to the 
gas kinetic temperature Tj^. The coupling coeffi- 
cient j/c oc rifi has been numerically evaluated by 
2]; it becomes strong when 1 + 5 > 5[(1 -I- z)/20]^ 
and is not important in the smooth IGM in the 



4 



relevant redshift range. The third term in equa- 
tion ||2Jl describes the Wouthuysen-Field effect, 
in which Lya pumping couples the spin temper- 
ature to the color temperature of the radiation 
field TLya t4 9 lA (which equals Tk in most astro- 
physical environments |16|). The physics of this 
mechanism have been described recently in |84j . 
Lya pumping effectively couples Ts and Tk when 
the radiation background at the Lya line center 
satisfies Jq. > 10~^^ erg cm~^ s~^ Hz^^ sr^^. In 
a low-density, neutral IGM, such a background 
would need to be present in order for the IGM to 
be visible in the 21 cm transition. 

3. A Short History of the Universe 

It is obvious from the above discussion that 
the IGM brightness temperature 5T{v) depends 
sensitively on the thermal history of the IGM. 
The spin temperature and optical depth there- 
fore depend on the kinetic temperature and den- 
sity of the IGM, as well as the mean intensity of 
the Lya flux. Once Thomson scattering of CMB 
photons becomes inefficient at the thermal decou- 
pling redshift Zd ^ 140, the IGM cools adiabati- 
cally until the first objects collapse. During this 
era, Tk < Tcmb; observations at < 30 MHz 
could in principle detect the IGM in absorption 
during this epoch \61\ . The cooling trend reverses 
itself as structure begins to form, but the sub- 
sequent temperature evolution is both inhomo- 
geneous and highly uncertain. While early es- 
timates suggested that Lya photons themselves 
would inject significant thermal energy into the 
IGM, jHj showed that this heating channel is in re- 
ality quite slow. Instead, X-rays (primarily from 
supernovae or accreting black holes), adiabatic 
compression, and structure formation shocks are 
likely to control the temperature evolution of the 
IGM. Because of this complex spin temperature 
evolution, 21 cm signatures can be divided into 
three separate eras. Here we will describe each of 
these phases in turn. The top panel of Figure 
shows a cartoon of the qualitative temperature 
evolution in a double reionization model. To fix 
ideas, we will refer to this Figure in the following 
discussion. We show two examples for Tk- The 
thin line includes only adiabatic cooling and pho- 



toionization heating (which we have applied at 
the maximal ionized fraction from each genera- 
tion for simplicity) . The thick line has additional 
heat sources before overlap. 



3.1. The First Structures 

The first phase occurs while luminous sources 
are rare; radiative feedback (including both heat- 
ing and Lya coupling) can therefore be ignored 
in most of the IGM. Although Tk < Icmb in 
this phase, we thus have Ts — Tcmb and most 
of the gas remains invisible. In Figure ^ this 
era corresponds to z > 25. However, note that 
only moderate overdensities are required in or- 
der to force Ts and Tcmb to depart from each 
other. One regime in which collisions are im- 
portant are "minihalos": collapsed objects that 
are not massive enough to cool and form stars. 
With mean overdensities S ~ 200, these objects 
are dense enough for collisions to quickly force 
the spin temperature to the virial temperature of 
the gas. They would therefore appear in emis- 
sion against the CMB. While an individual mini- 
halo (with M ^10^ M0) is much too small to 
be observed with the SKA, the integrated signal 
from the distribution of minihalos will cause rms 
fluctuations of order ST ~ 0.1-1 mK on angular 
scales of a few arcminutes ,27,, depending sensi- 
tively on z. 

However, as described in §2, collisional cou- 
pling does not require virial densities. More mod- 
erate overdensities, such as those that appear dur- 
ing the collapse of proto-sheets and fllaments, suf- 
flce to heat the gas and force Ts Tk- While 
these shocks are difficult to model without de- 
tailed numerical simulations, semi-analytic esti- 
mates suggest that shocked gas can amplify the 
signal expected from minihalos by a factor of a 
few at 2; ~ 20 J7|. Finally, bubbles of emis- 
sion or absorption may appear around the first lu- 
minous sources because their radiation fields can 
trigger the Wouthuysen-Field mechanism locally 
inSI- Thus, in this phase we expect relatively 
weak fluctuations to trace the formation of the 
flrst baryonic structures in the universe. 
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Figure 2. The brightness temperature of the 21 cm transition at several redshifts, as predicted by the 
"late reionization" simulation analyzed in |18| . Each panel corresponds to the same slice of the simulation 
box (with width 10h~^ comoving Mpc and depth = 0.1 MHz), at z = 12.1, 9.2, and 7.6, from left 
to right. The three epochs shown correspond to the early, middle, and late stages of reionization in this 
simulation; we assume Ts » Tqmb throughout. 



3.2. First Light 

This first phase ends when luminous sources 
become sufficiently numerous so as to make the 
Wouthuysen-Field effect important globally. At 
this point Ts Tk in the diffuse IGM and it 
becomes visible against the CMB Assuming 
that the Lya background turns on faster than the 
X-ray background, most of the IGM will still be 
cold and will initially be visible in absorption. In 
the cartoon of Figure ^ this era begins at z ~ 25 
as the spin temperature falls below Tcmb- Super- 
posed on the absorbing background, shock-heated 
filaments and minihalos will still appear in emis- 
sion. Thus in this epoch the 21 cm fluctuations 
depend on the density, temperature, and ionized 
fraction. The situation is sufficiently complex 
that predictions require detailed numerical sim- 
ulations (to trace the temperature structure from 
shocks in the IGM) as well as detailed model- 
ing of the radiation background (to describe the 
Lya coupling and X-ray heating). Such predic- 
tions are not yet available on the large (several 
Mpc) scales relevant for 21 cm observations, al- 



though simulations on smaller scales are now be- 
coming available Nevertheless, this phase 
is particularly interesting because it signals the 
appearance of the first luminous sources. For ex- 
ample, as shock-heating continues and the X-ray 
background builds, the diffuse IGM will be heated 
above Tcmb and the mean absorption will shift to 
emission, as occurs at z ~ 21 in Figure ^ Com- 
paring the onset of Lya coupling and the onset of 
heating will thus constrain the characteristics of 
the first luminous sources. 

3.3. A Uniformly Heated IGM 

Once the Lya and X-ray backgrounds become 
sufficiently strong, the 21 cm fluctuation pattern 
becomes independent of the temperature field. 
Even though the IGM will likely have a com- 
plex temperature field from continued shock heat- 
ing, the brightness temperature is essentially un- 
affected: T ^ 1 in equation The 21 cm signal 
therefore depends only on 6 and xh during this 
phase. As a result, this epoch is more straight- 
forward to consider from a theoretical standpoint 
than those described above. Several groups have 
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analyzed it recently with both numerical "9"18"21' 
and analytic [53iil9i20j techniques. In Figure ^ 
this approximation is accurate for all z < 20 or 
so. An example of the evolving fluctuation pat- 
tern during this phase is shown in Figure |2J 

The ultimate observational goal of the SKA is 
to map the evolving structure seen in the cen- 
tral panel of Figure [21 Structures in the IGM are 
causally related to the collapsed objects that are 
responsible for the photoionizing flux, and mea- 
suring these interelations will teach us about the 
first generations of stars and protogalaxies. How- 
ever, direct imaging of the IGM in the 21cm line 
is a formidable challenge, and it is likely that the 
first steps in observing the 21cm signal from the 
EOR will be statistical studies that discriminate 
between evolutionary scenarios through measure- 
ments of (for example) the angular power spec- 
trum. 

If uniform heating occurs sufficiently early com- 
pared to reionization, fluctuations in xh will also 
be negligible. In this case the 21 cm field traces 
the density and yields a direct measurement of the 
power spectrum of density fluctuations at z ~ 10 
|48) . Such an observation would provide a sensi- 
tive probe of the growth of fluctuations during the 
early universe. The scales accessible to observa- 
tions are linear at this time, so predictions of the 
fluctuation spectrum are relatively easy to make 
[48I5H| ■ Zaldarriaga et al. pointed out the math- 
ematical analogy between measurements of CMB 
fluctuations and the 21 cm pattern, with the key 
difference that (because it is a line transition), 21 
cm maps contain independent information at dif- 
ferent frequencies. They showed how to describe 
the 21 cm fluctuations in terms of their angular 
power spectrum, which can in turn be computed 
from the three-dimensional power spectrum -P(k). 
The fluctuations on a scale I — 2tt/9 are parame- 
terized by the Ci coefficients; the rms fluctuations 
are then given by ((5T2)1/2 ^ i{i + 1)Ci/2t: (note 
1 arcmin corresponds to Z « 2 x 10^). In the 
case in which density fluctuations dominate, this 
is simply the usual linear power spectrum on the 
relevant scales. The fluctuations increase with de- 
creasing angular scale, with a characteristic rms 
magnitude of several mK. The dotted curve in 
Figure |31 shows the angular power spectrum for 



this case of a universe with early, uniform heat- 
ing at z = 18 and with xh — 0.96. This spec- 
trum essentially registers the density fluctuations 
at z = 18 and has no strong features. 
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Figure 3. An example of the time evolution of the 
21 cm power spectrum. The time ordering is: dot- 
ted (z = 18, Xh = 0.96), short-dashed (z = 15, 
Xh = 0.81), long-dashed (z = 13, xh = 0.52) and 
solid (z = 12, Xh = 0.26). The thick curves as- 
sume the reionization model of jTH] while the thin 
dashed curve has each galaxy host an indepen- 
dent HII region (z = 13, xh — 0.52). The model 
assumes that each baryon inside a star forming 
galaxy produces ^ = 40 ionizing photons. From 



Once luminous sources become important, the 
primary source of fluctuations are the HII regions 
around ionizing sources, as shown in the middle 
panel of Figure El Predictions of the fluctuation 
pattern therefore depend on the particular model 
of reionization. One approach is to use numerical 
simulations to follow reionization. This is difficult 
for several reasons. First, radiative transfer re- 
mains a numerically challenging problem. It has 
usually only been added as a "post-processing" 
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step to existing hydrodynamic simulations. Sec- 
ond, each simulation represents a substantial in- 
vestment of computer time, and it is presently 
impractical to examine a large parameter space 
of reionization models. Third, locating typical 
ionizing sources (with M ^ 10^ M©) at these 
redshifts requires high mass resolution, limiting 
the size of the computational box to L < lOft.^^ 
comoving Mpc 0^ . This severely limits the abil- 
ity of numerical simulations to predict the 21 cm 
pattern on the arcminute scales (corresponding to 
several comoving Mpc) that would be accessible 
to the SKA. Nevertheless, much can be learned 
about the fluctuation pattern by examining simu- 
lations on small scales and extrapolating to larger 
scales |9I18I21| . These studies again suggest rms 
fluctuations 5T ~ 1-10 niK on scales of ~ 1-10 
arcminutes {I ^ lO'^-lO^). The fluctuations also 
depend on the spectral bandwidth of the obser- 
vation; they begin to be significantly suppressed 
when the radial scale (set by the spectral channel 
resolution) exceeds the transverse scale |53I19| . 
Thus the relevant angular scales correspond to 
- 0.2-2 MHz, or Az - 0.02-0.2. 
Analytic models present a different set of prob- 
lems. Most importantly, reionization is a highly 
nonlinear process, and constructing an accurate 
model for the sizes and spatial distribution of 
HII regions is challenging. The simulations de- 
scribed above all show that the usual approxima- 
tion - single HII regions surrounding each galaxy 
- does not describe overlap well. Instead, dense 
regions of the universe containing many galax- 
ies rapidly grow into large ionized regions (with 
scales of a few comoving Mpc) once the ionizing 
flux escapes from the over-densities into the low 
density regions. Thus analytic treatments must 
take into account large scale features of the den- 
sity field and not simply the regions near single 
galaxies. Moreover, there are a number of uncer- 
tain parameters (such as the dumpiness of gas in 
the IGM) that must be calibrated to simulations. 
Nevertheless, they offer the advantage of allowing 
efficient parameter studies and are not limited to 
small scales. Recently, JHl have developed an 
analytic model of reionization that includes the 
crucial features of the reionization process. The 
model relies on a simple extension of the Press- 



Schechter model 001 to compute the ionization 
pattern. It assumes that each galaxy can ionize 
a fixed multiple C, of its own baryonic mass. It 
then determines the sizes of HII regions by find- 
ing those regions that contain enough collapsed 
gas (i.e., galaxies) to ionize the entire region; [T^ 
showed that this can be approximated relatively 
simply through the "excursion set" formalism 
Given the background cosmology, the size distri- 
bution depends on only two parameters: the num- 
ber of ionizations per collapsed baryon {(,), and 
the minimum halo mass that can host an ionizing 
source. 

Given a model for reionization, we can then de- 
termine the observable characteristics of a 21 cm 
measurement. With high signal-to-noise maps, 
the size distribution of HII regions and the char- 
acteristics of reionization can be measured in 
a straightforward way. However, even without 
maps, statistical measurements are still extremely 
powerful. As described above, the simplest ap- 
proach is to compute the power spectrum of 21 cm 
emission, which depends on both the size distri- 
bution of HII regions and fluctuations in the back- 
ground density fleld. It is also relatively easy to 
interpret. Several examples of this simple statis- 
tic are shown in Figures 13 and 0] (both are taken 
from |19p. Figure shows how the power spec- 
trum evolves with time as the global neutral frac- 
tion decreases from unity to xh = 0.26. As the 
HII regions appear and grow, they add a feature 
to the power spectrum on scales of several ar- 
cminutes {I ~ lO^-lO*). The rms fluctuations in- 
crease by a factor of several on the characteristic 
scale of the HII regions, yielding a clear signa- 
ture of reionization even in statistical measure- 
ments. The crucial point is that the power spec- 
trum evolves rapidly through reionization, so its 
measurement will allow us to trace the time his- 
tory of the HII regions. The thin curve shows 
an example where the sizes of HII regions are de- 
termined by individual galaxies (as if the galaxies 
were randomly distributed); it emphasizes the im- 
portance of the model of reionization in predict- 
ing the power spectrum features and, conversely, 
the power of 21 cm tomography to differentiate 
these models. Note that these Figures assume 
infinite frequency resolution; as discussed above. 
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the fluctuations are significantly suppressed if tfie 
comoving radial distance exceeds the transverse 
distance |53I19| . Until late in reionization (when 
the HII regions become very large), the desired 
bandwidths correspond to Az ~ 0.02-0.2, and 
the redshift evolution within each bin is modest. 



10 - 




10 100 1000 lO" 10= 



Figure 4. The 21 cm power spectrum in three dif- 
ferent models of "double reionization." The solid 
curves assume that the universe was uniformly 
ionized to xh = 0.5 by a first generation of ob- 
jects. The others assume that the first generation 
shut off while discrete bubbles still existed, again 
when Xh = 0.5. The second generation of sources 
completes reionization in all cases. In each model, 
the thin curves have a total xh = 0.46 and the 
thick curves have xr = 0.27. From ^0) . 



Figure 01 shows the kind of detailed informa- 
tion that we can extract from 21 cm tomography 
through the power spectrum |20j . Here we con- 
trast predictions for two different prescriptions of 
"double reionization," a scenario first considered 
in j6l5()| . Specifically, we consider the mecha- 
nism through which the first generation of highly 
ionizing sources shuts off. One possibility (solid 
curves) is that these sources ionize the universe 



more or less uniformly (or that they fully ionize 
the universe, which then recombines to a nearly 
uniform neutral fraction). Another possibility is 
that the first generation shuts off while its HII 
regions still exist, so that the neutral fraction is 
highly inhomogeneous during the transition. The 
two sets of dashed curves in Figure 01 show this 
prescription, with two different sets of properties 
for the first generation of sources. We see that, 
in the former case, the power spectrum and the 
bubble feature are significantly suppressed. In 
the latter scenario, however, the bubble feature 
is quite strong. Essentially, the first generation 
imprints a characteristic bubble size that persists 
until the second generation of sources are able to 
dominate the total number of ionizations. While 
these scenarios are no doubt overly simplified, 
they serve to illustrate the power of 21 cm ob- 
servations. 

It is important to note, however, that the power 
spectrum is only the simplest statistic available, 
essentially measuring the rms deviation as a func- 
tion of angular scale. There are several other ways 
to approach the problem that take expicit advan- 
tage of the three-dimensional nature of the ob- 
servations |37] . A related question is how to take 
advantage of the non-gaussian nature of the sig- 
nal. Because fluctuations in the linear density 
field are gaussian, the power spectrum provides 
a complete characterization of that field. How- 
ever, the ionized fraction is most assuredly not a 
gaussian random field; it is actually closer to a 
discrete variable (taking values zero or unity) in 
the simplest models. During the middle to late 
stages of reionization, other statistical properties 
may be better descriptors of the observations. For 
example, argue that such statistics can help 
to distinguish the way in which reionization pro- 
ceeds, whether from high-to-low density regions 
or from low-to-high density regions. More de- 
tailed methods to extract this information from 
the measurements are still required, however. 
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4. Observational Challenges & Instrument 
Requirements 

4.1. Foreground Contamination 

The surface brightness fluctuations 5T ^ 1- 
10 mK in the redshifted 21cni spectral hue 
amount to ^^10"^ of the brightness of the radio 
sky, which is dominated by the diffuse contin- 
uum synchrotron emission from the Milky Way 
Galaxy. The redshift range 6.2 < z < 20 corre- 
sponds to l%7 > v> 68 MHz [where v ^ 1420(1+ 
z)~^ MHz]. In this frequency range, a good rule 
of thumb for the brightness temperature in the 
coldest directions is ~ 180(i^/180MHz)-2-SK. 
The greatest complication imposed by the Galac- 
tic foreground is the increase in receiving sys- 
tem noise Tgys, since, according to the radiometer 
equation, the rms noise in the antenna temper- 
ature AT ~ Tgys/V^^ tint rises in proportion 
to sky brightness {Tsys > Tsky); achieving a low 
noise level must be accomplished by increasing 
integration time Unt, since the bandwidth Aiy is 
set by matching the expected widths of the spec- 
tral features (~0.2 to 2 MHz). The brightness 
temperature of the radio sky increases toward 
the Galactic plane and is an order of magnitude 
higher toward the Galactic Genter, necessitating 
the identification of selected high Galactic lati- 
tude regions in the sky for studying 21 cm emis- 
sion from the EOR. The Galactic synchrotron 
halo is considered to be smooth on arcminute 
scales when mapped at normal sensitivities. How- 
ever, the halo's fluctuations at mK sensitivity re- 
main to be explored. 

To flrst order, conventional radio spectroscopic 
techniques for continuum subtraction should be 
capable of removing the foregrounds created by 
the Galactic synchrotron and by discrete extra- 
galactic sources J5J. Interferometric antenna ar- 
rays are insensitive to smooth emission, and "re- 
solve away" the DC term of the Fourier trans- 
form of the sky brightness. This leaves behind the 
structure caused by fluctuations in the Galactic 
foreground and the discrete radio source popula- 
tions on the sky. Gontinuum source spectra are 
typically either power laws or vary only slowly 
from a power law. As a consequence, there are 
two steps to radio continuum subtraction, de- 



pending on the strength of the contaminating 
sources. 

In the flrst step, the coordinates and spectrum 
of the strongest sources are measured, and the in- 
strumental response to these sources is subtracted 
directly from the calibrated u-v dataset. This 
leaves behind the many weaker sources, whose 
sum provides the continuum emission that pop- 
ulates the 1 to 20 arcminute resolution pixels in 
which we wish to measure the 21cm EOR emis- 
sion. 




K.A, 



Figure 5. Spectral line data cube. The spectral 
information for one R.A, Declination coordinate 
has been drawn in for emphasis. Since radio in- 
terferometers do not measure the 'zero spacing 
flux density, ' the average flux density of each fre- 
quency plane of an observation cube is equal to 
zero. Spectral baseline removal involves fitting 
the residual continuum for each R.A-Dec pixel 
separately with a low-order polynomial. The fit- 
ted function is then subtracted from the values as 
a function of frequency, (from ^) 



The second step operates on a pixel by pixel 
basis. The standard continuum subtraction oper- 
ates on the intermediate stage data-product from 
spectral line observations, which is a data-cube 
(see Figure 0) composed of many image planes, 
with each plane of different frequency. The sums 
of the power law spectra with different spectral in- 
dices that fill a R.A. -Declination pixel may have 
spectral curvature, but these are adequately fit 
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by second order polynomials. Several authors 
[Tl>3 8 10 12 have raised concerns about the con- 
fusing effects of discrete source foregrounds, but 
conventional techniques are adequate for their re- 
moval, as illustrated in Figure 

Other approaches to the problem of continuum 
removal have also been studied. They too rely 
on the spectral smoothness of the contaminants. 
One way to approach it is through the symmetries 
in the data cube in Fourier space PZ] . Another is 
through frequency differencing |53j . where neigh- 
boring planes in the data cube are subtracted 
from each other. This technique is similar to 
that described above except that the baseline sub- 
traction is performed implicitly on the visibilities. 
The top panel of Figure |7| illustrates the level to 
which the foregrounds can be removed in the for- 
malism of |53j . The solid and dotted curves show 
estimates of the 21 cm fluctuation spectrum at 
z — 10, with and without HII regions. The dot- 
dashed curve shows the level to which the radio 
point source population predicted by can be 
removed. It is clear that smooth foregrounds al- 
low for more than the necessary precision in con- 
tinuum removal. 

Additional uncertainties surround the sub- 
tleties of instrumental response to the known (and 
yet to be discovered) foregrounds. For example, 
an array telescope actually senses every contin- 
uum source above the horizon at some level; the 
response has strong frequency dependence that 
increases with angular separation of each source 
from the telescope pointing coordinates. The dis- 
tribution of sources (in strength and location) 
produces a pseudo-random superposition of these 
effects in the net response and could create spuri- 
ous spectral features that will confuse 21cm line 
studies of the EOR. The solution may be to per- 
form the data reduction for EOR surveys with ac- 
cess to all-sky radio source catalogues including 
coordinates, flux densities and spectral shapes. 
Another subtlety arises from the unknown polar- 
ization properties of the Galactic halo foreground, 
in which differential Faraday rotation of the syn- 
chrotron emission will enter in a frequency depen- 
dent way. If the telescope polarization response 
is not calibrated to high precision, then the rotat- 
ing polarization vector could couple in the analy- 
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Figure 6. An illustration of spectral baseline fore- 
ground subtraction for 3'x3' pixels. Top panel: 
average integral flux density vs. frequency per 
pixel; this includes radio point sources but as- 
sumes that both the smooth Galactic contribu- 
tion and the continuum sources brighter than 
Scut = 100 /iJy have been subtracted. Bottom 
panel: Open triangles are rms deviations with- 
out continuum subtraction. Open squares are 
rms residuals after linear spectral baseline sub- 
traction. Filled circles are rms residuals after 
quadratic baseline subtraction. In each case, the 
upper curve includes fluctuations from large scale 
structure , while the lower curve includes only 
Poisson statistical fluctuations (based on mean 
source count statistics). The horizontal dotted 
line marks the 10 mK level at 150 MHz, which is 
typical of the rms fluctuations expected for the 
EoR signal. The righthand axis indicates logTe 
for this pixel size (from [3]). 
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sis to the measured intensity, and this could vary 
with frequency and position in the sky to mimic 
EOR signals. Much remains to be done in the 
exploration of instrumental effects, astrophysical 
foregrounds, and their interactions within the ob- 
servational process. 



1 1 

\ 0.1 



!' 1 1 


^ ': 




^ ' 
























1 1 1 



10 100 103 



Figure 7. Instrument sensitivities, as estimated 
by 1231 • In both panels the dotted and solid 
curves show simple estimates for the 21 cm sig- 
nal with and without reionization, respectively. 
The top panel shows the noise power spectrum 
(relevant for imaging) for LOFAR (long-dashed 
curve) and SKA (short-dashed curve). The bot- 
tom panel shows estimates for the errors in loga- 
rithmic /-bins of a statistical measurement of the 
power spectrum. The dot-dashed curve in the top 
panel shows the estimated level of residual fore- 
ground contamination from faint radio continuum 
sources. 



4.2. SKA requirements 

The goal of observing 21 cm emission from the 
EOR constrains the frequency coverage and array 
configuration of the telescope. For example, the 
SKA must be able to observe at frequencies on 
the order of 100 MHz. The opacity in the Gunn- 
Peterson trough in high redshift quasars suggests 
that reionization may be ending at z ~ 6 [1^ . If 
the universe is largely neutral (xu > 0.3) up to 
this epoch, even limiting the frequency coverage 
to > 150 MHz will provide valuable information 
about reionization. On the other hand, the CMB 
data presently indicates that reionization must 
begin at z ~ 15 [2H|- If the universe is partially 
ionized beyond this point, the SKA would need to 
push to much lower frequencies in order to make 
21 cm measurements. Unfortunately, precise con- 
straints on the reionization history will probably 
not appear for several years (see §5) and may even 
be best determined by 21cm line observations. In 
any case, it appears that measurements of the 
formation of the first structures, as described in 
§3.1-3.2, could require frequency coverage below 
70 MHz. 

Sensitivity to faint fluctuations in surface 
brightness against the bright radio sky at these 
frequencies will require a large collecting area. 
The achievable uncertainty in surface brightness 
ATb w ATa/r]f = Tsys/rif^/ Av Unt for emission 
that fills the array beam depends on the array 
filling factor ijf. Because the fluctuation spec- 
trum falls off toward large angular scales, the ar- 
ray must have diameter of ^3000 wavelengths if 
the beam is to match 1' structure. On the other 
hand, to obtain surface brightness sensitivity, the 
aperture cannot be too severely diluted due to the 
rjf dependence. 

The filling factor requirement, along with the 
frequency dependence of the radio sky, are il- 
lustrated in Figure |S1 This Figure provides an 
overview of the issues that influence the ability of 
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Figure 8. Estimated detection limits of faint EOR features at four redshifts. The brightness fluctuation 
spectrum for a ACDM model assuming pure density fluctuations (following US|) has been scaled upward 
by a factor of three to provide a rough measure of the 3-cr peaks of brightness. We assume lS.v = 0.5 MHz 
and neglect redshift space distortions. Solid lines indicate 5a detection limits for 100 and 1000 hour SKA 
observations, for different diameter arrays obtained by diluting 1 km^ of collecting area to achieve higher 
angular resolution. Dotted curves do the same for a LOFAR array. 
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the SKA to map specific features in the EOR, 
which wiU be crucial for relating structures in 
the IGM to star forming regions observed at in- 
frared and possibly millimeter wavelengths. The 
Figure has four panels in order to show the in- 
creasing degree of difficulty of the observation 
toward higher redshifts: as the frequency de- 
creases, the sky brightness and the diffraction 
limited resolution worsen. The signal predictions 
assume pure density ffuctuations (i.e., they ne- 
glect the contribution from HII regions) and that 
Ts >> TcMB everywhere. They are generated 
as described in JH! (who in turn followed j48j ) . 
assuming Af = 0.5 MHz and neglecting redshift 
space distortions (which will amplify the signal 
by a small factor). Note that we show Scr peaks 
in order to illustrate the brightest features. The 
instrumental sensitivity is best for the largest an- 
gular scales (> 10') since the array beam would 
be filled by the brightness fluctuations for a filled 
array (r]f — 1). For finer angular scales, the aper- 
ture must be diluted in order to match the beam 
to the fluctuations, causing the sensitivity to de- 
teriorate. For this figure, the noise level is set by 
the Tsys ~ Tsky , with the assumption that residu- 
als from foreground subtraction make a negligible 
contribution. 

The statistical measures of the fiuctuation 
power spectrum have somewhat different depen- 
dencies. The short-dashed curve in the upper 
panel of Figure [7| shows the approximate noise 
power spectrum for the SKA assuming one month 
of continuous observing, a 0.4 MHz channel, and 
Acs/Tsys = 2 X lO^m^/ K (see reference for 
details). This curve essentially shows the noise 
per visibility as a function of scale and is analo- 
gous to the noise per pixel in a real-space map. 
The long-dashed curve shows the approximate 
noise for the Low Frequency Array (LOFAR). We 
can see clearly that (because the sky brightness 
is essentially fixed) high signal-to-noise maps re- 
quire collecting areas on the order of a square 
kilometer; even then they will be difficult. We 
note that this sensitivity curve was constructed 
with several fairly crude approximations and is 
only accurate to a factor of a few. It depends in 
detail on the configuration of the array, the an- 
tenna beams, the correlator properties, and the 



observing strategy. In the future, these questions 
must be addressed in more detail. Also note that, 
on the scales of the HII regions, the contrast be- 
tween ionized and neutral regions will be substan- 
tially larger than the simple rms variation shown 
here; in this sense map-making during the middle 
stages of reionization will be easier than suggested 
by this plot. 

Of course, statistical properties are much eas- 
ier to measure than maps, because we can average 
over many measurements on each scale. The re- 
quirements for the signal-to-noise ratio on each 
angular scale are therefore much weaker. The 
bottom panel of Figuredshows the estimated er- 
rors in such a measurement for the SKA and LO- 
FAR, assuming logarithmic bins in Z-space and a 
100 square degree field of view. We see that the 
power spectrum of the fluctuations can be mea- 
sured with good accuracy even with collecting ar- 
eas smaller than a square kilometer provided that 
the field of view is large. 

Finally, the most important angular scales for 
measuring the fluctuations are in the range ~ 1- 
20', which contain both the peak of fluctuations 
in the linear density fleld and the expected scale 
of HII regions during most of reionization (at least 
in the model of ^Hl)- This suggests that the SKA 
have a compact core with baselines 3 km that 
contains a substantial fraction of the collecting 
area (Figure H assumes that 20% of the antennae 
are inside a core of diameter 3 km). While long 
baselines will be useful for point source removal, 
21 cm tomography requires heavy emphasis on 
the compact core. 

5. Conclusions 

Understanding reionization is one of the major 
goals of cosmology. We have argued that 21 cm 
tomography is a particularly powerful way to at- 
tack the problem, but there are several others. 
We are already reaching the limits of observa- 
tions with the Lya forest; observations of com- 
plete Gunn-Peterson troughs at z ~ 6 indicate 
that this tool cannot be taken to substantially 
higher redshifts. Moreover, Lya absorption is so 
strong that even measurements of zero transmis- 
sion do not imply a fully neutral medium. Thus, 
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more detailed constraints from optical observa- 
tions must rely on modeling of the HII region 
around quasars [51135) . measurements of the de- 
tailed shape of the Lya absorption profile 36' in 
(for example) gamma-ray burst afterglow spectra, 
or the careful interpretation of Lya emission lines 
from high-redshift galaxies 139132141122171 . While 
invaluable, these methods are fraught with uncer- 
tainty. 

More detailed measurements of the CMB will 
also yield more information about the ionization 
history. Most importantly, WMAP and Planck 
will better measure the total optical depth to 
electron scattering and place a strong integral 
constraint on the ionization history. However, 
the CMB power spectrum is relatively insensi- 
tive to the evolution of the ionized fraction |25j . 
and strong constraints will probably have to await 
space-borne missions beyond Planck. (Note that 
the reionization signal appears on larger angular 
scales than ground-based or balloon experiments 
can easily probe.) 

Finally, direct observations of high-redshift 
protogalaxies, with (for example) the James 
Webb Space Telescope or the upcoming gener- 
ation of twenty-meter ground-based telescopes, 
will obviously teach us an enormous amount 
about the growth of the first luminous objects. 
However, we will still have only indirect measure- 
ments of the IGM structure and the process of 
reionization. 

Thus it seems that 21 cm tomography pro- 
vides a unique and invaluable probe of the high- 
redshift universe. No other technique can mea- 
sure the three-dimensional distribution and time 
evolution of neutral gas in such detail. Also, only 
21 cm measurements can extend into the era be- 
fore the formation of the first luminous objects 
(see §3.1 and 3.2) to map the initial stages of 
structure formation. Provided that the technical 
challenges described in §4 can be overcome, there 
is no doubt that 21 cm tomography will revolu- 
tionize our understanding of reionization and the 
early universe. While statistical measurements 
of the distribution of ionized gas may be possi- 
ble with LOFAR or some other instrument be- 
fore SKA, we expect that a collecting area near 
one square kilometer will be required in order to 



make high signal-to-noise maps of the neutral gas 
distribution. 

On the other hand, we expect that these mea- 
surements will complement the next generation 
of infrared. X-ray, and CMB instruments. For 
example, comparing the distribution of ionized 
gas to the luminous sources in a cosmological vol- 
ume would show explicitly how those sources are 
responsible for ionization. We can also use 21 
cm observations of large ionized regions to locate 
bright quasars and in principle to constrain many 
of their unknown properties, including their ages 
and spectra [tHj . 

Finally, we note that reionization is an evolv- 
ing field. The data are continuing to improve, and 
our theoretical understanding is advancing even 
more rapidly. At the same time, the amount of 
attention focused on 21 cm studies has expanded 
dramatically, and our understanding of the in- 
strumental and data analysis challenges is im- 
proving by leaps and bounds. This review has 
presented a snapshot of the field in May 2005. 
Our expectations will no doubt change in many 
ways (some unexpected) before the SKA begins 
observations. 
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